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Abstract 

We introduce a new approach to the maximum flow problem in undirected, capacitated graphs using a- 
*K^ congestton-approximators: easy-to-compute functions that approximate the congestion required to route 

^_^ single-commodity demands in a graph to within a factor of a. Our algorithm maintains an arbitrary 

f^ flow that may have some residual excess and deficits, while taking steps to minimize a potential function 

^S| measuring the congestion of the current flow plus an over-estimate of the congestion required to route 

» I the residual demand. Since the residual term over-estimates, the descent process gradually moves the 

^ I contribution to our potential function from the residual term to the congestion term, eventually achieving 

a flow routing the desired demands with nearly minimal congestion after 0{ae~^ log^ n) iterations. Our 
approach is similar in spirit to that used by Spielman and Teng (STOC 2004) for solving Laplacian 
systems, and we summarize our approach as trying to do for £oo-flows what they do for ^2-flows. 

Together with a nearly linear time construction of a n°'^'-congestion-approximator, we obtain 1 -f 

r/^ £-optimal single-commodity flows undirected graphs in time m}'^°^^' £~^ , yielding the fastest known 

algorithm for that problem. Our requirements of a congestion-approximator are quite low, suggesting 

even faster and simpler algorithms for certain classes of graphs. For example, an a-competitive oblivious 

Y\ routing tree meets our definition, even without knowing how to route the tree back in the graph. For 

I I graphs of conductance 4>, a trivial ^"^-congestion-approximator gives an extremely simple algorithm for 

finding 1 + e-optimal- flows in time 0{rn(p~^). 
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The maxiinum flow problem and its dual, the minimum cut problem are fundamental combinatorial opti- 
mization problems with a wide variety of applications. In the well-known maximum s — t flow problem we 
^T are given a graph G with edge capacities Ce, and aim to route as much flow as possible from s to i while 

restricting the magnitude of the flow on each edge to its capacity. We will prefer instead to think in terms of 
the equivalent problem of routing a single unit of flow from s to t while minimizing the maximum congestion 
L" l/e/ce| on any edge; clearly the minimum congestion for unit flow is equal to one divided by the maximum 

. !^ flow of congestion one. Once formulated that way, we need no longer restrict ourselves to s — t flows; given 

S^ a demand vector b € M" specifying the excess desired at each vertex, we aim to find a flow / e M™ with 

Jh divergence equal to b that minimizes the maximum congestion |/e/ce|r| 

In this paper, we introduce a new approach to this problem in undirected graphs. We maintain a flow 
that may not quite route b exactly, but we also keep track of an upper bound on how much it will cost 
us in congestion to fix it back up. We will aim to minimize a potential function measuring the current 
congestion plus an over-estimate on the cost of fixing up the residuals. By not needing to worry about 
precisely conserving flow at every vertex, we can take large steps in each iteration towards minimizing our 
potential function. On the other hand, by intentionally over-estimating the cost of fixing up the residuals, 
in the course of minimizing our potential function we must inevitably fix them up, as it will cost strictly less 
to do so. 

For a graph G, let C be the m x m diagonal matrix containing the edge capacities, and let B be the 
n X m divergence matrix, where {Bf)i is the excess at vertex i. For a set S* C ]/, we'll write bs — X^ies ^5'' 
the total excess in S, and cs = '^e-s^v\s '^ei the capacity of the cut (S*, V\S) in G. A valid demand vector 
satisfies by = 0- 
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^In fact the s — t case is no less general, since one could always add a new vertices s and t, connect each ti to s or t according 
to the sign of b„ with an edge of capacity /3|6u| and scale /9 until the additional edges are saturated. 



The minimum congestion flow problem for demands 6, and its dual, the maximum congested cut, are 

min IIC-i/lloo s.t.Bf = b (1) 

max b'^v si. ||CB^i;||i < 1 (2) 

We refer to the optimum value of these problems as opt (6). It is well-known that for problem (l2|, one of the 
threshold cuts with respect to v achieves bs/cs > b^v. 

1.1 History 

Much of the early work on this problem considers the general, directed edge case, culminating in the still- 
best binary blocking flow algorithm of Goldberg and Rao[l] that achieves 0(mmin(m^/^, n"^/^) time. Karger 
and Levine^ give evidence that the undirected case seems easier in graphs with small flow values. The 
smooth sparsiflcation technique of Benczur and Karger [1] shows one can split a graph with m edges into 
t — 0{me'^/n) graphs each with with 0{ne~^) edges, and each of which has at least (1 — e)/t of the capacity 
of the original graph. Therefore, for undirected graphs, any algorithm running in time T{m, n) can be 
replaced with one running in time 0{m) + (m/n)T{On, n). Using the algorithm of Goldberg and Rao yields 
approximate maximum flows in 0{m-n}/'^) time, and for many years that was the best known. 

In a breakthrough, Christiano et al. show how to compute approximate maximum flows in 0{mn^i^) 
time[3]. Their new approach uses the nearly linear-time Laplacian solver of Spielman and Teng[12, to 
take steps in the minimization of a softmax approximation of the edge congestions. Each step involves 
minimizing ||VFC~^/||2, a weighted ^2-norm of the congestions. While a naive analysis yields immediately 
yields a method that makes 0{\/rn) such iterations (because || • II2 approximates || • ||oo by a \/rn factor), 
they present a surprising and insightful analysis showing in fact only 0{'m}/'^) such (.2 iterations suffice. The 
maximum-flow-speciflc parts of ^3J are quite simple, needing only to maintain the weights W and then using 
the Spielman- Teng solver as a black box. 

1.2 Outline 

In this work we pry open that box, extract the parts we need, and apply them directly to the maximum 
flow problem. In the course of doing so, we push the running time for this decades-old problem nearly 
down to linear. Solving a Laplacian system Lv* = 6*, where L = BCB^ , corresponds to finding a fiow /* 
with BJ* = b* minimizing ||C~^"/||2. While Spielman and Teng work entirely in the dual space of vertex 
potentials and never explicitly represent flows, we can still translate between spaces throughout the algorithm 
to get an idea of what is actually going on. At any point, the vertex potentials v induce an optimal flow for 
some demands b — Lv; just perhaps not the ones desired. The solution of Spielman and Teng is to maintain 
a simpler graph G" that approximates G, in the sense that the £2 cost of routing flows in G" is within some 
factor a of the cost in G. The residual flow is routed optimally in G" (recursively), by solving L'v' — b* — b. 
The potentials that induced that flow are added to v, in the hope that Lv + Lv' ~ Lv + L'v' = b* . Indeed, 
if G" approximates G by a factor a, then b* — b gets smaller in a certain norm by (1 — l/a), nudging the 
flow towards actually routing b* . 

Our first step towards obtaining a ST-like algorithm is the definition of a good approximator for the 
congestion required by ^00-fiows. 

Definition 1.1. An a -congestion- approximator for G is a matrix R such that for any demand vector b, 

\\Rh\\oo < opt(6) < a||i?&||oo 

Our main result is that we can use good congestion approximators to quickly find near-optimal flows in 
a graph. We prove the following theorem in section [2j 

Theorem 1.2. There is an algorithm that, given demands b and access to an a -congestion- approximator 
R, makes 0{ae^^log (n)) iterations and then returns a flow f and cut S with Bf = b and ||G~"'^/|| < 
(1 -|- e)bs/cs. Each iteration requires 0{m) time, plus a multiplication by R and R^ . 

For the sake of giving the reader a concrete example of what a congestion-approximator might look like 
before continuing, we'll begin with two simple toy examples. 



Example 1.3. Let T he a maximum weight spanning tree in G, and let R be the n — 1 x n matrix with a 
row for each edge in T , and 

{Rb)e = - 

where (S*, V \ S) is the cut in G induced by removing e from T. 
Then, R is a m-congestion-approximator. 

Proof. Since {Rb)e is the congestion on the cut in G induced by removing e from T, certainly opt{b) > ||-Rt'||oo- 
On the other hand, at least 1/ni of the capacity of those cuts is contained in T, so routing b through T 
congests e by at most m\{Rb)e\- Multiplication by R and R^ can be done in 0{n) time via elimination on 
leaves. D 

For graphs of large conductance, we can obtain a trivial approximator by simply looking at how much 
the demand into each vertex congests its total degree. 

Example 1.4. Let G have conductance <j). Let R be n x n diagonal matrix with Ri^i — l/deg(z) where 
deg(«) — C{ij.. Then, R is a cj)^^ -congestion-approximator. 

Proof. Routing \bi\ into or out of vertex i certainly must congest one of its edges by at least |&j|/deg(i), so 
opt(6) > |ji?6||oo- On the other hand, the capacity of any cut in G is at least (j) times the total degree of the 
smaller side. It follows that if no vertex is congested by more than /3, then no cut is congested by more than 

Those two simple examples are analogous to the simple cases for £2 flows in the ST-algorithm. The 
former is analogous to preconditioning by a small-stretch spanning tree, while the latter is analogous to not 
preconditioning at all. As in the ST-algorithm, those two simple examples in fact capture the ideas behind 
our real constructions. In section pj we show how to apply the j-tree decomposition of MadrvflO] (itself 
based on the ultrasparsifiers of Spielman and Teng[T2])) to obtain good congestion approximators for any 
graph. 

Theorem 1.5. For any 1 < fc < logn, there is an algorithm to construct a data structure representing a 
log{n)^''' -congestion- approximator R in time 0{m -\- n^^^''^). Once constructed, R and R^ can each be 
applied in time 0{n^'^^'''). 



Combining theorem 1.5 where k = 6{\/logn) with theorem 1.2 yields our main result of a nearly-linear 
time algorithm for minimum congestion flows. 

Theorem 1.6. There is an algorithm to compute (1 -f- e)- approximate minimum congestion flows in time 
me~'^ ■ exp(0(yTogri)). For graphs of conductance (j), the same can be done in time 0{m(f>^^). 

2 Congestion Potential 

The key to our scheme lies transforming problem ([I]) to an unconstrained optimization problem, using the 
congestion-approximator to bound the cost of routing the residual. To that end, we introduce our potential 
function. 

min||C-Vl|oo + 2a||i?(6-S/)||oo (3) 

We believe the mere statement of the potential function (|3| to be the most important idea in this paper. 
Once ([3]) has been written down, the remainder of our algorithm is nearly obvious. We solve problem ([3]) 
nearly-optimally by approximating || • ||oo with a softmax. The softmax is well-approximated by its gradient in 
the region where steps are taken with ^oo-norm 0(1). Since ||-Ri?C||oo-foo ^ 1 for a congestion-approximator, 
we can take steps on edge e of size ^{a~^)ce. We include the details at the end of this section, proving the 
following theorem. 

Theorem 2.1. There is an algorithm AlmostRoute that given b and s < 1/2, performs 0{ae^'^logn) 
iterations and returns a flow f and cut S with, 

IIC-Vlloo + 2a\\R{b - B/)|U <{l + e) — 



Each iteration requires 0{m) time plus a multiplication by R and R 
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The flow / may not quite route the demands we wanted. Fortunately, that will be easy to fix. The extra 
factor of two in equation ([3| means that half of the contribution to the objective value from the residual 
part is pure slack. On the other hand, if / is nearly optimal, there can't be too much slack. 

Lemma 2.2. Suppose ||C-i/l|oo + 2a\\R{b - Bf)\\^ < (1 + e)opt{b). Then, \\R{b - Bf)\\^ < £||i?6||oo. 

Proof. Let / meet the assumption. Let /' be a routing of 6 — Bf in G with ||C^^/'||oo < Q;||i?(6 — i3/)||oo- 
Then, moving from / to / + /' decreases the objective value by atleast a\\R{b — Bf)\\. On the other hand, 
that decrease can't exceed eopt(6). Since opt(6) < a||i?6||oo, the lemma follows. D 

So while AlmostRoute may not route b, our bound for the congestion to route the residual is at most 
half of our bound to route the original demands. Furthermore, the objective already pays the cost of routing 
that residual. In fact, it pays it with a factor of two, so we need only route the remaining residual within a 
factor-two of optimal. That suggests an obvious way to route demands b: repeatedly invoke AlmostRoute 
on the remaining residual, until the congestion required to route it is extremely small compared to the 
congestion required to route b. Then route the final residual in a naive way, such as via a maximal spanning 
tree. The cost of that final routing will be paid for by the slack in the objective value of the first routing, 
simply by finding a factor 3/2-optimal routing for each residual after the first case. 



Formalizing the latter argument completes our proof of theorem 1.2 Set ba -s— b, and let {fo,So) ^— 
AlmostRoute (feo,£)- Next for i = 1,...,T where T — log(2TO), set bi <— 6i_i — B/i-i and {fi,Si) <— 
AlmostRoute (fei, 1/2) (we don't actually need any Si after So). Finally, let br+i — h — Bft, and let /t+i be 
a flow routing br+i in a maximal spanning tree of G. Output /i + • • ■ + fr+i and 5*0. Observe that theorem 
[2T] yields 

\\C-^fo\\oo + 2a\\Rb4^ < (l+e)bsJcso 
IIC-V.IU + 2a||i?&,+i||oo < (3/2)opt(6,0 < (3/2)a|li?6,|U 

Beginning with the former inequality and repeatedly applying the latter yields, 

(1 + e)^ > \\C-'M\oo + 2a\\Rbi\\^ 

> (l/2)a||i?6i||oo + IIC^Volloo + ||C7-Vi||oo + 2a||i?62||oo 

> (l/2)a||i?6i||oo + \\C-'fo\\oo + ■■■ + ||C-Vt||oo + 2a||i?&T+i||oo 
On the other hand, by choice of T, we have 

||C"Vt+i||oo < ma||i?6T+i||oo < ma2-^|li?5i|loo < {l/2)a\\Rbi\\^ 
Combining the two yields the theorem. 

2.1 Proof of Theorem \2A\ 



For this preliminary draft, we prove theorem 2.1 with slightly worse parameters of 0(a^e~'^log (n)), using 
naive steepest descent. The better parameters follow from using Nesterov's accelerated gradient method [TT] 
instead, as will be proved in the final version of this paper. Of course, for the case of a = exp(0(\/logn), 
then still a^ = exp(0(-v/logn)), so this naive analysis still yields nearly linear-time flow algorithms. 
We approximate || • ||oo using the symmetric softmax function. 



lmax(a;) = log ^ e"^' -I- e "^^ ] 



We make use of some elementary facts about Imax. 
Fact 2.3. Letx,y^W^. Then, 



||Vlmax(a;)||i < 1 (4) 

Vlmax(a;)^a; > lmax(a;) — log(2(i) (5) 

|Vlmax(a;) - Vlmax(2/)||i < \\x - y\\oo (6) 



We will approximate problem ([3| with the potential function, 

0(/) = lmax(C" V) + Imax {2aR{b - Bf)) (7) 

Since equation (l7| approximates equation ([3]) to within an additive 0{\ogn), we will be concerned with 
minimizing 0(/) after scaling /, b so 0(/) = 9{e~^ logn). 



AlmostRoute(&, e): 

• Initialize f — 0, scale b so 2Q;||i?6||oo = 16£^^ log(?i). 

• Repeat: 

- While (j){f) < We-^ log(n), scale / and b up by 17/16. 

- Set,5^||CV0(/)||i. 

- If (5 > e/4, set fe^fe- i+i^ sgn(V</.(/)e)ce 

- Otherwise, terminate and output / together with the potentials induced by V(/)(/)(see below), 
after undoing any scaling. 



Each step requires computing V0(/), which requires 0{m) time plus a multiplication by R and a mul- 
tiplication by R^ . Further, the partial derivative of the residual part for a particular edge is equal to a 
potential difference between the cndpoints of that edge. When (f>{f) is nearly-optimal, those potentials yield 
a good dual solution for our original problem. 

Lemma 2.4. When AlmostRoute terminates, we have a flow f and potentials v with, 

IIC-i/lloo + 2a\\R{b - B/)|U < (1 + e) ^ " 



\CB^v\\i 

Proof. Set xi = C^^f, X2 — 2aR{b~Bf), andp^ — Vlmax(a;i). Set v ~ R^P2 to be our potentials. Observe 
that V(/)(/) = C^^pi - 2aB^v. First, equation Q yields 

2a\\CB'^v\\i < ||pi||i + ||pi-2aCB^^;||i <l + d 

By equation ^, using the fact that C and R have at most n^/2 rows and </)(/) > 16e^^ logn, 

pJC-\f + 2aplR{b - Bf) > <P{f) - Alogn > 0(/)(l ~ e/4) 

On the other hand, 

s<f>{f) > ||cv<^(/)||i||c-Vlloo 

= {C-'pi-2aB^R^P2ff 
= pJC-^f~2apjRBf 

Combining the two yields, 

2ab'^v> (/>(/)(l-£/4-(5) 

Altogether, using the fact that 5 < e/4 at termination, we have 

b'^v ^ 0(/)(l-£/2) ^ (/>(/) 

\\CB^v\\i - l + e/4 - 1 + e 

Observing that </>(/) overestimates ||C^"'^/||oo + 2a||i?(6 — -B/)||oo completes the proof. D 

Lemma 2.5. AlmostRoute terminates after at most 0(a^e~"^ logn) iterations. 

Proof. Let us call the iterations between each scaling a phase. Since ||i?6||oo gives us the correct scale to 
within factor a, we will scale at most O(loga) times. 



Lethe — — Y:rfjj2 sgn(V/0(/)e)ce be our step. Then, equation (Ig]), together with the fact that ||_Ri3C||oo-i-oo < 
1 for a congestion-approxiniator i? yields, 

</-(/ + M < 0(/)+V0(/r/i+i^^||C-i/i||L 

= 0(/) - r!(e2«-2) 

Since we raised 0(/) by at most e^^ logn when scaling, and each step drops i^(/) by at least r2(e^a^^), 
there can be at most 0{o?e~'^\o%n) steps between phases. D 

3 Computing Congestion-Approximators 



In this section we prove theorem 1.5 using a construction of MadrvfTO]. itself based on a construction of 
Spielman and Teng[12 . 

Definition 3.1 (Madry lOj). A j-tree is a graph formed by the union of a forest with j components, together 
with a graph H on j vertices, one from each component. The graph H is called the core. 

Theorem 3.2 (Madrv[TU]'l. For any graph G and t > 1, we can find in time 0{tm) a distribution oft graphs 
(Ai, Gi) such that, 

• Each Gi is a O (m log m/t) -tree, with a core containing at most m edges. 

• Gi dominates G on all cuts. 

• ^j AiGi can be routed in G with congestion O(logn). 



We briefly remark that while the statement of theorem 3.2 in IQJ contains an additional logarithmic 
dependence on the capacity-ratio of G, that dependence is easily eliminated. We elaborate further in appendix 
\K\ Our construction will simply apply theorem 3.2 recursively, sparsifying the core on each iteration. To 



accomplish that, we use an algorithm of Bencziir and Karger[l . 

Theorem 3.3 (Bencziir, Karger Ij). There is an algorithm Sparsify(G, e) that, given a graph G with m 
edges, takes 0{m) time and returns a graph G' with in' — 0{ne^^ logn) edges such that the capacity of cuts 
in the respective graphs satisfy 

G<G' < (1 + £)G 

Further, the edges of G' are scaled versions of a subset of edges in G, with no edge scaled by more than 
(1 + £)m/m' . 

We now present the algorithm for computing the data structure representing a congestion-approxiniator. 
The algorithm ComputeTrees assumes its input is sparse; our top-level data-structure is constructed by 
invoking ComputeTrees(Sparsif y(G, 1), n^'^), where k is the parameter of theorem 



1.5 



ComputeTrees(G, t): 


• If n = 1, return. 


• Using theorem 


3.2 


compute distribution (Ai,Gi)* of max(l, n/i)-trees. 


• Pick the t graphs of largest A^, throw away the rest, and scale the kept A^ to sum to 1. 


• For i = l,...,t: 


— H[ -i— Sparsify(i/i, 1), where Hi is the core of Gi 


— Li i- ComputeTrees(iJj',i) 


• Return the list L — (A.^, Fi, Li)j^i where Fi is the forest of G,;. 



The analysis of ComputeTrees correctness will make use of another algorithm for sampling trees. The 
SampleTree procedure is only used for analysis, and is not part of our flow algorithm. 



SampleTree(L = (A^, F^, Ti)*): 

• Pick i with probability A^. 

• Output Fi + SampleTree (Ti) 



Lemma 3.4. Let G have 0{n) edges, and set L -^ ComputeTrees(G, t). Then, every tree in the sample 
space 0/ SampleTree (L) dominates G on all cuts, and E[SampleTree(_L)] is routable in G with congestion 
log(n)'°s(")/'o9(*). Further, the computation of L takes 0{tn) time. 

Proof. By induction on n. For n — 1 the claim is vacuous, so suppose n — t'^'^^. Since G has O(nlogn) edges, 
the distribution output by theorem 3.2 will have 0(t log^ n) entries. We have H'^ > Hi and Hi + Fi > G. 



Furthermore, the inductive hypothesis implies that every tree Ti in ScmipleTree(Li) dominates H^. Then, 

T, + F,> H' + F, > H, +F, = G.,>G 



Sparsifying the distribution from 0(ilog n) to t scales Xi by at most 0(log n), so that X]i=i -^j^j i^ 
routable in G with congestion at most log n larger than the original distribution. Since H^ < 2Hi, by 
the multicommodity max-flow/min-cut theorem "91 H'i is routable in Hi with congestion O(logri). By the 
inductive hypothesis, E[SampleTree(Li)] is routable in Hi with congestion log ' (n). It follows then that 
E[SampleTree(L)] is routable in G with congestion at most log ' ^^'(n). 

Finally, ComputeTrees requires 0(tn) time to compute the distribution, another 0{tn) time to sparsify 
the cores, and then makes t recursive calls on sparse graphs with ji/t vertices. It follows that the running 
time of ComputeTrees is 0{tn). D 

Lemma 3.5. Let R be the matrix that has a row for each forest edge in our data structure, and {Rb)e is the 
congestion on that edge when routing b. //E[SampleTree(i)] is routable in G with congestion a, then R is 
a a-congestion-approximator for G. Further, R has 0{tn) rows. 

Proof. Since the capacity of each tree-edge dominates the capacity of the corresponding cut in G, opt(6) > 
||i?6||oo- On the other hand, b can be routed in every tree with congestion ||i?6||oo- By routing a Pr[T] 
fraction of the flow through tree T, we route b in E[SampleTree(_L)] with congestion ||i?6||oo- But then b can 
be routed in G while congesting by at most an a factor larger. 

The total number of edges in R satisfies the recurrence E{n) < nt + tE(n/t) as each edge is either in one 
of the t toplcvcl forests, or in one of the t subgraphs. D 

Having constructed our representation of R, it remains only to show how to multiply by R and R^ . We 
use the following lemmas as subroutines, which are simple applications of leaf-elimination on trees. 

Lemma 3.6. There is an algorithm TreeFlow that, given a tree T and a demand vector b, takes 0{n) time 
and outputs for each tree edge, the flow along that edge when routing b in T. 

Lemma 3.7. There is an algorithm TreePotential that, given a tree T annotated with a price Pe for each 
edge, takes 0{n) time and outputs a vector of vertex potentials v such that, for any i,j, the sum of the prices 
on the path from j to i in T is Vi — Vj . 

We begin with computing R. We take as input the demand vector b, and then annotate each forest edge 
e with the congestion r,, induced by routing b through a tree containing e. 



Coinputei?(5, i = (A„F„T,)*): 
• For i = 1, . . . ,t: 

— Let T be the tree formed by taking Fi, adding a new vertex s, and an edge from s to each 
core- vertex of Fi. Augment b with demand zero to the new vertex. 

— /^ TreeFlow(6,T). 

— Set Te <— fe/ce for cach forest edge in Fi. 

— Set 6' to a vector indexed by core- vertices, with 6' equal to the flow on the edge from s to 
core- vertex j. 

— ComputeR(6',Tj). 



Lemma 3.8. The procedure Compute R{b, L) correctly annotates each edge e with r^ ~ {Rb)e, and takes 
0{tn) time. 

Proof. Let L = {\i,Fi,Ti)\. We argue by induction on the depth of recursion. Fix a level and index i. 
Observe that the cut in G induced by cutting a forest edge is the same regardless of what tree T lies on 
the core: it is the cut that separates the part of Fi not containing the core from the rest of the vertices. It 
follows that we may place any tree on the core vertices, invoke TreeFlow, and obtain the flow on each forest 
edge. Next, for each component S oi Fi, the total excess bs must enter 5* via the core vertex. It follows that 
in a flow routing b on Fi + T' , for any tree T', the restriction of that flow to T' must have excess bs on the 
core vertex of S, so it suffices to find a flow in the core with demands 6'- = bsj ■ But routing 6 in i^^ + T will 
place exactly 65. units of flow on the edge from s to core- vertex j. 

The running time consists of t invocations of TreeFlow cach taking 0{n) time, plus t recursive calls on 
graphs of size n/t, for a total running time of 0{tn). D 

To compute i?^, we assume each forest edge e has been annotated with a price Pe that must be paid by 
any flow per unit of congestion on that edge, and output potentials v such that Vi — Vj is the total price to 
be paid for routing a unit of flow from j to i. 



CompnteR^ iL = {X,,F,,T,)l): 

• w <— 

• For i — 1, . . . ,t: 

— v' ^ Computei?^(ri). 

— Let T be the tree formed by taking Fi, adding a vertex s, and an edge from s to each core- vertex 
of Fi. Set q^ = Pe/C(. for each forest edge, and q^ = f' for edge e from s to core- vertex j . 

— v" <s— TreePotential(T, g) 

— Add v" to V after removing the entry for s. 

• Return v 



Lemma 3.9. Given edges annotated with per- congestion prices, the procedure Computei?^(L) correctly re- 
turns potentials v .such that Vk — Vj is the cost per unit of flow from vertex j to k. 

Proof. Let L = [Xi, Fi,Ti)\. We argue by induction on the depth of recursion. Fix a level; a flow must pay 
its toll to each Gi, so the resulting potential equals the sum of the potentials for each i. Fix an index i. A 
unit of flow from j to k is first routed from j to the core- vertex of the component of Fi containing j, then to 
the core-vertex of the component containing k, and then finally to k. By induction, we assume that v' yields 
potentials that give the per-unit costs of routing between core-vertices. Placing a star on the core with the 
edge from s to core-vertex j having per-unit cost w' preserves those costs. If p^ is the price of an edge per 
unit of congestion, then q^. = Pe/c^ is the price of an edge per unit of flow. It follows that the total toll paid 
is the same as the toll paid in T; thus, the potentials output by TreePotential(r, q) are correct. 

The running time consists of t recursive calls to Computei?^ on graphs of size n/t, plus t invocations of 
TreePotential each taking 0{n) time, for a total running time of 0(tn). D 



4 Final Remarks 

We remark that there are many other ways to obtain good congestion approximators. The obhvious routing 
schemes of [SI [2] require polynomial time to compute, but, once computed, give us a single tree whose single- 
edge cuts yield a log(n)'^^^)-congestion approximator. Furthermore, we only need the actual tree, and not 
the routings of the tree back in the original graph. If such a single tree could be computed in nearly-linear 
time, it would make an ideal candidate for use in our algorithm. 

There have been substantial simplifications to Spielman and Teng's original algorithm (see 0[3)- It may 
be possible to use some of those techniques to further simplify our algorithm. 
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A Fixing Theorem 3.2 



The proof of theorem 3.2 maintains a length function /(e) for each edge, and repeatedly invokes an algorithm 



SmallStretchTree that returns a spanning tree T on G with, 

E^T(e)c(e)<0(logn)^/(e)c(e) (8) 



where Irie) is the length of the path between e's endpoints in the tree. Without loss of generality, by scaling, 
we assume '^f,l{e)c{e) — m. Let xl^ic') = 1 if e' is a tree edge that lies on the path in T containing e. 
Then, 

J2 h{e)c{e) = ^ c(e) ^ cM{e, e'W) = ^ l{e') ^ x(e, e')c{e) = ^ l{e')cT{e') 

e e e' e' e e' 

where crie') is the total capacity of edges routed through e' in T. 

The dependence on the capacity ratio arises from the fact that there may be many different scales of 
congestion on the edges of T. The solution is simply to replace / with ^'(e) — ^ — ; ^ mixture of the 

original lengths with the inverse capacities. Constructing a small-stretch tree with respect to I' still satisfies 
equation (|8| with an extra factor of two, but also implies no tree edge is congested by more than 0{m). 
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